Phylogenomic Studies Research Articles

Discovery of a new Pseudalomya Telenga, 1930 (Hymenoptera, Ichneumonidae, Ichneumoninae) species from Taiwan and its implications for the systematic position of this genus

The rare genus Pseudalomya Telenga comprises two species, which are found only in the Eastern Palaearctic region and high mountains of the Oriental region. The phylogenetic position of Pseudalomya remains unclear because of its intermediate morphology between two ichneumonine tribes, Alomyini and Phaeogenini. This article reports the discovery of a new species of Pseudalomya: Pseudalomya truncaticornissp. nov. Specimens were collected during a survey of insect fauna in the Dasyueshan area of Shei-Pa National Park, one of the high-altitude regions in Taiwan. The new species can be diagnosed by its body coloration, frontal horn shape, facial punctures, metasomal tergite sculpture, and wing venation. To the best of our knowledge, this is the first record of Pseudalomya in Taiwan. This article also presents a diagnostic key to the global species of Pseudalomya. In this study, molecular phylogenetic analyses were performed using one mitochondrial and two nuclear gene sequences from P. truncaticornissp. nov. and other members of the Ichneumoniformes group. The results indicate that Pseudalomya should be classified within Phaeogenini, distinct from Alomyini, but more comprehensive phylogenomic studies are needed to confirm this placement.

Kaistella rhinocerotis sp. nov., isolated from the faeces of rhinoceros and reclassification of Chryseobacterium faecale as Kaistella faecalis comb. nov.

Strain Ran72T, a novel Gram-stain-negative, obligately aerobic, non-motile, and rod-shaped bacterium, was isolated from the faeces of the rhinoceros species Ceratotherium simum. The novel bacterial strain grew optimally in Reasoner's 2A medium under the following conditions: 0 % (w/v) NaCl, pH 7.5, and 30 °C. Based on phylogenetic analysis using 16S rRNA gene sequencing, strain Ran72T was found to be most closely related to Chryseobacterium faecale F4T (98.4 %), Kaistella soli DKR-2T (98.0 %), and Kaistella haifensis H38T (97.4 %). A comprehensive genome-level comparison between strain Ran72T with C. faecale F4T, K. soli DKR-2T, and K. haifensis H38T revealed average nucleotide identity, digital DNA-DNA hybridization, and average amino acid identity values of ≤74.9, ≤19.3, and ≤78.7 %, respectively. The major fatty acids were anteiso-C15 : 0 (22.3 %), with MK-6 being the predominant respiratory quinone. The major polar lipids of strain Ran72T were phosphatidylethanolamine, four unidentified aminolipids, and two unidentified lipids. Based on our chemotaxonomic, genotypic, and phenotype characterizations, strain Ran72T was identified as representing a novel species in the genus Kaistella, for which the name Kaistella rhinocerotis sp. nov. is proposed, with the type strain Ran72T (=KACC 23136T=JCM 36038T). Based on the outcomes of our phylogenomic study, Chryseobacterium faecale should be reclassified under the genus Kaistella as Kaistella faecalis comb. nov.

Characterization of the complete chloroplast genome sequence of <i>Hypericum petiolulatum</i> (Hypericaceae)

<i>Hypericum</i> L. (Hypericaceae) is one of the best-selling herbal medicines in the world comprising ca. 500 species of herbs, shrubs, and small trees. <i>Hypericum petiolulatum</i> Hook. f. & Thomson ex Dyer is widely distributed in China, Vietnam, Myanmar, Nepal, India, Malaysia, and Bhutan and is used as a traditional herb to treat hemoptysis and inflammation. In this study, we sequenced and assembled the complete chloroplast (cp) genome of <i>H. petiolulatum</i>. The complete plastome of <i>H. petiolulatum</i> was 136,105 bp in length, with a large single copy region (LSC) of 93,709 bp, a small single copy region (SSC) of 11088 bp and two identical inverted repeats (IRs) of 15,654 bp. The overall GC content of the plastome was 37.0%, while GC contents of the LSC, SSC, and each IR were 35.5%, 31.0%, and 43.8%, respectively. In addition, 116 genes consisting of 76 protein-coding genes, six ribosomal RNA genes, and 34 transfer RNA genes were identified. A phylogenetic analysis of 14 taxa inferred based on cp genome sequences revealed a close relationship between <i>H. petiolulatum</i> and <i>H. perforatum</i>. The complete cp genome sequence of <i>H. petiolulatum</i> reported in this paper will facilitate population and phylogenomics studies of this medicinal plant group.

MAST: Phylogenetic Inference with Mixtures Across Sites and Trees.

Hundreds or thousands of loci are now routinely used in modern phylogenomic studies. Concatenation approaches to tree inference assume that there is a single topology for the entire dataset, but different loci may have different evolutionary histories due to incomplete lineage sorting, introgression, and/or horizontal gene transfer; even single loci may not be treelike due to recombination. To overcome this shortcoming, we introduce an implementation of a multi-tree mixture model that we call MAST. This model extends a prior implementation by Boussau et al. (2009) by allowing users to estimate the weight of each of a set of pre-specified bifurcating trees in a single alignment. The MAST model allows each tree to have its own weight, topology, branch lengths, substitution model, nucleotide or amino acid frequencies, and model of rate heterogeneity across sites. We implemented the MAST model in a maximum-likelihood framework in the popular phylogenetic software, IQ-TREE. Simulations show that we can accurately recover the true model parameters, including branch lengths and tree weights for a given set of tree topologies, under a wide range of biologically realistic scenarios. We also show that we can use standard statistical inference approaches to reject a single-tree model when data are simulated under multiple trees (and vice versa). We applied the MAST model to multiple primate datasets and found that it can recover the signal of incomplete lineage sorting in the Great Apes, as well as the asymmetry in minor trees caused by introgression among several macaque species. When applied to a dataset of four Platyrrhine species for which standard concatenated maximum likelihood and gene tree approaches disagree, we observe that MAST gives the highest weight (i.e. the largest proportion of sites) to the tree also supported by gene tree approaches. These results suggest that the MAST model is able to analyse a concatenated alignment using maximum likelihood, while avoiding some of the biases that come with assuming there is only a single tree. We discuss how the MAST model can be extended in the future.

Robust demarcation of the family Peptostreptococcaceae and its main genera based on phylogenomic studies and taxon-specific molecular markers.

The family Peptostreptococcaceae, which contains 15 genera including Clostridioides, presently lacks proper circumscription. Using 52 available genomes for Peptostreptococcaceae species, we report comprehensive phylogenomic and comparative analyses to reliably discern their evolutionary relationships. In phylogenetic trees based on core genome proteins and 16S rRNA gene sequences, the examined species formed a strongly supported clade designated as Peptostreptococcaceae sensu stricto. This clade encompassed the genera Peptostreptococcus (type genus), Asaccharospora, Clostridioides, Intestinibacter, Paeniclostridium, Paraclostridium, Peptacetobacter, Romboutsia and Terrisporobacter, and two misclassified species (viz. Eubacterium tenue and 'Clostridium dakarense'). The distinctness of this clade is strongly supported by eight identified conserved signature indels (CSIs), which are specific for the species from this clade. Based on the robust evidence provided by presented studies, we are proposing the emendment of family Peptostreptococcaceae to only the genera within the Peptostreptococcaceae sensu stricto clade. We also report 67 other novel CSIs, which reliably demarcate different Peptostreptococcaceae species clades and clarify the classification of some misclassified species. Based on the consistent evidence obtained from different presented studies, we are making the following proposals to clarify the classification of Peptostreptococcaceae species: (i) transfer of Eubacterium tenue, Paeniclostridium ghonii and Paeniclostridium sordellii as comb. nov. into the genus Paraclostridium; (ii) transfer of Clostridioides mangenotii as a comb. nov. into Metaclostridioides gen. nov.; (iii) classification of 'Clostridium dakarense' as a novel species Faecalimicrobium dakarense gen. nov., sp. nov. (type strain FF1T; genome and 16S rRNA accession numbers GCA_000499525.1 and KC517358, respectively); (iv) transfer of two misclassified species, Clostridium paradoxum and Clostridium thermoalcaliphilum, into Alkalithermobacter gen. nov.; and (v) proposals for two novel families, Peptoclostridiaceae fam. nov. and Tepidibacteraceae fam. nov., to accommodate remaining unclassified Peptostreptococcaceae genera. The described CSIs specific for different families and genera provide novel and reliable means for the identification, diagnostics and biochemical studies on these bacteria.

Challenging old microbiological treasures for natural compound biosynthesis capacity.

Strain collections are a treasure chest of numerous valuable and taxonomically validated bioresources. The Leibniz Institute DSMZ is one of the largest and most diverse microbial strain collections worldwide, with a long tradition of actinomycetes research. Actinomycetes, especially the genus Streptomyces, are renowned as prolific producers of antibiotics and many other bioactive natural products. In light of this, five Streptomyces strains, DSM 40971T, DSM 40484T, DSM 40713T, DSM 40976T, and DSM 40907T, which had been deposited a long time ago without comprehensive characterization, were the subject of polyphasic taxonomic studies and genome mining for natural compounds based on in vitro and in silico analyses. Phenotypic, genetic, and phylogenomic studies distinguished the strains from their closely related neighbors. The digital DNA-DNA hybridization and average nucleotide identity values between the five strains and their close, validly named species were below the threshold of 70% and 95%-96%, respectively, determined for prokaryotic species demarcation. Therefore, the five strains merit being considered as novel Streptomyces species, for which the names Streptomyces kutzneri sp. nov., Streptomyces stackebrandtii sp. nov., Streptomyces zähneri sp. nov., Streptomyces winkii sp. nov., and Streptomyces kroppenstedtii sp. nov. are proposed. Bioinformatics analysis of the genome sequences of the five strains revealed their genetic potential for the production of secondary metabolites, which helped identify the natural compounds cinerubin B from strain DSM 40484T and the phosphonate antibiotic phosphonoalamide from strain DSM 40907T and highlighted strain DSM 40976T as a candidate for regulator-guided gene cluster activation due to the abundance of numerous "Streptomyces antibiotic regulatory protein" (SARP) genes.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution.

Angiosperms (flowering plants) are by far the most diverse land plant group with over 300,000 species. The sudden appearance of diverse angiosperms in the fossil record was referred to by Darwin as the "abominable mystery," hence contributing to the heightened interest in angiosperm evolution. Angiosperms display wide ranges of morphological, physiological, and ecological characters, some of which have probably influenced their species richness. The evolutionary analyses of these characteristics help to address questions of angiosperm diversification and require well resolved phylogeny. Following the great successes of phylogenetic analyses using plastid sequences, dozens to thousands of nuclear genes from next-generation sequencing have been used in angiosperm phylogenomic analyses, providing well resolved phylogenies and new insights into the evolution of angiosperms. In this review we focus on recent nuclear phylogenomic analyses of large angiosperm clades, orders, families, and subdivisions of some families and provide a summarized Nuclear Phylogenetic Tree of Angiosperm Families. The newly established nuclear phylogenetic relationships are highlighted and compared with previous phylogenetic results. The sequenced genomes of Amborella, Nymphaea, Chloranthus, Ceratophyllum, and species of monocots, Magnoliids, and basal eudicots, have facilitated the phylogenomics of relationships among five major angiosperms clades. All but one of the 64 angiosperm orders were included in nuclear phylogenomics with well resolved relationships except the placements of several orders. Most families have been included with robust and highly supported placements, especially for relationships within several large and important orders and families. Additionally, we examine the divergence time estimation and biogeographic analyses of angiosperm on the basis of the nuclear phylogenomic frameworks and discuss the differences compared with previous analyses. Furthermore, we discuss the implications of nuclear phylogenomic analyses on ancestral reconstruction of morphological, physiological, and ecological characters of angiosperm groups, limitations of current nuclear phylogenomic studies, and the taxa that require future attention.

Actinoplanes oblitus sp. nov., producing the glycopeptide antibiotic A477.

In 1973, Eli Lilly and Company described the filamentous actinomycete producing the glycopeptide antibiotic A477 as an Actinoplanes species on the basis of its morphological and physiological features and deposited it as NRRL 3884T. In this paper, we report that the phylogenetic analysis based on the 16S rRNA gene sequence and the whole genome phylogenomic study indicate that NRRL 3884T forms a distinct monophyletic line within the genus Actinoplanes, being most closely related to Actinoplanes octamycinicus NBRC 14524T [99.6 % 16S rRNA gene similarity, 89.4 % average nucleotide identity (ANI), 46.0 % digital DNA-DNA hybridization (dDDH)] and Actinoplanes ianthinogenes NBRC 13996T (98.8 % 16S rRNA gene similarity, 89.0 % ANI, 47.0 % dDDH). NRRL 3884T forms an extensively branched, non-fragmented vegetative mycelium; either sterile aerial hyphae or regular subglobose sporangia are formed depending on cultivation conditions. The cell wall contains meso-2,6-diaminopimelic acid and 2,6-diamino-3-hydroxypimelic acid and the diagnostic sugars are glucose, mannose and ribose with a minor amount of rhamnose. The predominant menaquinone (MK) is MK-9(H4), with minor amounts of MK-9(H2), MK-9(H6) and MK-9(H8). Mycolic acids are absent. The diagnostic phospholipids are diphosphatidylglycerol and phosphatidylethanolamine. The major cellular fatty acids are anteiso-C17 : 0, iso-C16 : 0 and iso-C15 : 0, with moderate amounts of anteiso-C15 : 0 and iso-C17 : 0. The genomic G+C content is 71.5 mol%. Significant differences in the genomic, morphological, chemotaxonomic and biochemical data between NRRL 3884T and the two most closely related Actinoplanes type strains clearly demonstrate that NRRL 3884T represents a novel species of the genus Actinoplanes, for which the name Actinoplanes oblitus sp. nov. is proposed. The type strain is NRRL 3884T (=DSM 116196T).

Lessons from assembling UCEs: A comparison of common methods and the case of Clavinomia (Halictidae).

Sequence data assembly is a foundational step in high-throughput sequencing, with untold consequences for downstream analyses. Despite this, few studies have interrogated the many methods for assembling phylogenomic UCE data for their comparative efficacy, or for how outputs may be impacted. We study this by comparing the most commonly used assembly methods for UCEs in the under-studied bee lineage Nomiinae and a representative sampling of relatives. Data for 63 UCE-only and 75 mixed taxa were assembled with five methods, including ABySS, HybPiper, SPAdes, Trinity and Velvet, and then benchmarked for their relative performance in terms of locus capture parameters and phylogenetic reconstruction. Unexpectedly, Trinity and Velvet trailed the other methods in terms of locus capture and DNA matrix density, whereas SPAdes performed favourably in most assessed metrics. In comparison with SPAdes, the guided-assembly approach HybPiper generally recovered the highest quality loci but in lower numbers. Based on our results, we formally move Clavinomia to Dieunomiini and render Epinomia once more a subgenus of Dieunomia. We strongly advise that future studies more closely examine the influence of assembly approach on their results, or, minimally, use better-performing assembly methods such as SPAdes or HybPiper. In this way, we can move forward with phylogenomic studies in a more standardized, comparable manner.

Phylogenomics studies and molecular markers reliably demarcate genus Pseudomonas sensu stricto and twelve other Pseudomonadaceae species clades representing novel and emended genera.

Genus Pseudomonas is a large assemblage of diverse microorganisms, not sharing a common evolutionary history. To clarify their evolutionary relationships and classification, we have conducted comprehensive phylogenomic and comparative analyses on 388 Pseudomonadaceae genomes. In phylogenomic trees, Pseudomonas species formed 12 main clusters, apart from the "Aeruginosa clade" containing its type species, P. aeruginosa. In parallel, our detailed analyses on protein sequences from Pseudomonadaceae genomes have identified 98 novel conserved signature indels (CSIs), which are uniquely shared by the species from different observed clades/groups. Six CSIs, which are exclusively shared by species from the "Aeruginosa clade," provide reliable demarcation of this clade corresponding to the genus Pseudomonas sensu stricto in molecular terms. The remaining 92 identified CSIs are specific for nine other Pseudomonas species clades and the genera Azomonas and Azotobacter which branch in between them. The identified CSIs provide strong independent evidence of the genetic cohesiveness of these species clades and offer reliable means for their demarcation/circumscription. Based on the robust phylogenetic and molecular evidence presented here supporting the distinctness of the observed Pseudomonas species clades, we are proposing the transfer of species from the following clades into the indicated novel genera: Alcaligenes clade - Aquipseudomonas gen. nov.; Fluvialis clade - Caenipseudomonas gen. nov.; Linyingensis clade - Geopseudomonas gen. nov.; Oleovorans clade - Ectopseudomonas gen. nov.; Resinovorans clade - Metapseudomonas gen. nov.; Straminea clade - Phytopseudomonas gen. nov.; and Thermotolerans clade - Zestomonas gen. nov. In addition, descriptions of the genera Azomonas, Azotobacter, Chryseomonas, Serpens, and Stutzerimonas are emended to include information for the CSIs specific for them. The results presented here should aid in the development of a more reliable classification scheme for Pseudomonas species.

Compositae-ParaLoss-1272: A complementary sunflower-specific probe set reduces paralogs in phylogenomic analyses of complex systems.

A family-specific probe set for sunflowers, Compositae-1061, enables family-wide phylogenomic studies and investigations at lower taxonomic levels, but may lack resolution at genus to species levels, especially in groups complicated by polyploidy and hybridization. We developed a Hyb-Seq probe set, Compositae-ParaLoss-1272, that targets orthologous loci in Asteraceae. We tested its efficiency across the family by simulating target enrichment sequencing in silico. Additionally, we tested its effectiveness at lower taxonomic levels in the historically complex genus Packera. We performed Hyb-Seq with Compositae-ParaLoss-1272 for 19 Packera taxa that were previously studied using Compositae-1061. The resulting sequences from each probe set, plus a combination of both, were used to generate phylogenies, compare topologies, and assess node support. We report that Compositae-ParaLoss-1272 captured loci across all tested Asteraceae members, had less gene tree discordance, and retained longer loci than Compositae-1061. Most notably, Compositae-ParaLoss-1272 recovered substantially fewer paralogous sequences than Compositae-1061, with only ~5% of the recovered loci reporting as paralogous, compared to ~59% with Compositae-1061. Given the complexity of plant evolutionary histories, assigning orthology for phylogenomic analyses will continue to be challenging. However, we anticipate Compositae-ParaLoss-1272 will provide improved resolution and utility for studies of complex groups and lower taxonomic levels in the sunflower family.

Pyricularia Are Mostly Host-Specialized with Limited Reciprocal Cross-Infection Between Wheat and Endemic Grasses in Minas Gerais, Brazil.

Wheat blast, caused by Pyricularia oryzae Triticum (PoT), is an emerging threat to global wheat production. The current understanding of the population biology of the pathogen and epidemiology of the disease has been based on phylogenomic studies that compared the wheat blast pathogen with isolates collected from grasses that were invasive to Brazilian wheat fields. In this study, we performed a comprehensive sampling of blast lesions in wheat crops and endemic grasses found in and away from wheat fields in Minas Gerais. A total of 1,368 diseased samples were collected (976 leaves of wheat and grasses and 392 wheat heads), which yielded a working collection of 564 Pyricularia isolates. We show that, contrary to earlier implications, PoT was rarely found on endemic grasses, and, conversely, members of grass-adapted lineages were rarely found on wheat. Instead, most lineages were host-specialized, with constituent isolates usually grouping according to their host of origin. With regard to the dominant role proposed for signalgrass in wheat blast epidemiology, we found only one PoT member in 67 isolates collected from signalgrass grown away from wheat fields and only three members of Urochloa-adapted lineages among hundreds of isolates from wheat. Cross-inoculation assays on wheat and a signalgrass used in pastures (U. brizantha) suggested that the limited cross-infection observed in the field may be due to innate compatibility differences. Whether or not the observed level of cross-infection would be sufficient to provide an inoculum reservoir, or serve as a bridge between wheat growing regions, is questionable and, therefore, deserves further investigation.

An insight on the complete chloroplast genome of Gomphocarpus siniacus and Duvalia velutina, Asclepiadoideae (Apocynaceae).

We studied the complete chloroplast genome of Gomphocarpus siniacus and Duvalia velutina from Asclepiadoideae subfamily; due to their medicinal importance and distribution worldwide their interest became high. In this study we analyzed the complete chloroplast genomes of G. siniacus and D. velutina using Illumina sequencing technology. The sequences were compared with the other species from Apocynaceae family. The complete genome of G. siniacus is 162,570 bp while D. velutina has154, 478 bp in length. Both genomes consist of 119 genes; encode 31 tRNA genes, and eight rRNA genes. Comparative studies of the two genomes showed variations in SSR markers in which G. siniacus possesses 223 while D. velutina has 186. This could be used for barcoding in order to aid in easy identification of the species. Phylogenetic analysis on the other hand reaffirms the tribal position of G. siniacus in Asclepiadeae and D. velutina in Ceropegieae. These findings could be used in subsequent research studies of angiosperms identification, genetic engineering, herb genomics and phylogenomic studies of Apocynaceae family.

A comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important species.

Haplotype-resolved chromosomal-level genome assembly of Buzhaye (Microcos paniculata)

Microcos paniculata is a shrub used traditionally as folk medicine and to make herbal teas. Previous research into this species has mainly focused on its chemical composition and medicinal value. However, the lack of a reference genome limits the study of the molecular mechanisms of active compounds in this species. Here, we assembled a haplotype-resolved chromosome-level genome of M. paniculata based on PacBio HiFi and Hi-C data. The assembly contains two haploid genomes with sizes 399.43 Mb and 393.10 Mb, with contig N50 lengths of 43.44 Mb and 30.17 Mb, respectively. About 99.93% of the assembled sequences could be anchored to 18 pseudo-chromosomes. Additionally, a total of 482 Mb repeat sequences were identified, accounting for 60.76% of the genome. A total of 49,439 protein-coding genes were identified, of which 48,979 (99%) were functionally annotated. This haplotype-resolved chromosome-level assembly and annotation of M. paniculata will serve as a valuable resource for investigating the biosynthesis and genetic basis of active compounds in this species, as well as advancing evolutionary phylogenomic studies in Malvales.

Thoracic and abdominal outgrowths in early pterygotes: a clue to the common ancestor of winged insects?

One of the fundamental questions in insect evolution is the origin of their wings and primary function of ancestral wing precursors. Recent phylogenomic and comparative morphological studies broadly support a terrestrial ancestor of pterygotes, but an aquatic or semiaquatic ancestor cannot be ruled out. Here new features of the branchial system of palaeodictyopteran larvae of several different instars of Katosaxoniapteron brauneri gen. et sp. nov. (Eugereonoidea) from the late Carboniferous collected at Piesberg (Germany) are described, which consist of delicate dorsolateral and lamellate caudal abdominal gills that support an aquatic or at least semiaquatic lifestyle for these insects. Moreover, the similar form and surface microstructures on the lateral abdominal outgrowths and thoracic wing pads indicate that paired serial outgrowths on segments of both tagmata presumably functioned as ancestral type of gills resembling a protopterygote model. This is consistent with the hypothesis that the wing sheaths of later stage damselfly larvae in hypoxic conditions have a respiratory role similar to abdominal tracheal gills. Hence, the primary function and driving force for the evolution of the precursors of wing pads and their abdominal homologues could be respiration.

Kakapo: easy extraction and annotation of genes from raw RNA-seq reads.

kakapo (kākāpō) is a Python-based pipeline that allows users to extract and assemble one or more specified genes or gene families. It flexibly uses original RNA-seq read or GenBank SRA accession inputs without performing global assembly of entire transcriptomes or metatranscriptomes. The pipeline identifies open reading frames in the assembled gene transcripts and annotates them. It optionally filters raw reads for ribosomal, plastid, and mitochondrial reads, or reads belonging to non-target organisms (e.g., viral, bacterial, human). kakapo can be employed for targeted assembly, to extract arbitrary loci, such as those commonly used for phylogenetic inference in systematics or candidate genes and gene families in phylogenomic and metagenomic studies. We provide example applications and discuss how its use can offset the declining value of GenBank's single-gene databases and help assemble datasets for a variety of phylogenetic analyses.

Phylogenomic analysis of Bupleurum in Western Sichuan, China, including an overlooked new species

A comparative analysis of chloroplast (cp) genomes and 45s nuclear ribosomal DNA (nrDNA), and a phylogenomic study of six closely related species (including an overlooked new species) of genus Bupleurum from the western part of Sichuan Province in southwestern China were performed. The six species are similar morphologically and it is difficult to identify them; moreover, their genetic relationships remain unclear. It was found that the cp genomes of the six Bupleurum species were extremely similar, and they were highly homogeneous in terms of cp genome structure, genes and its arrangement. Intergenic spacer rpl32-trnL, petA-psbJ, trnK-rps16, and the coding gene ycf1 were considered highly variable. In phylogenetic trees constructed based on the complete cp genome, protein-coding sequences, nrDNA and ITS sequences, Chinese Bupleurum species all formed two major clades; among these trees, nrDNA tree had the best species resolution; the highly variable regions showed no advantage over other molecular markers. Among the six Bupleurum species, B. malconense, B. sichuanense were close relatives to B. chinense and B. yinchowense, B. chaishoui may also be a consanguinity, while B. microcephalum, B. wenchuanense, and the new species B. pseudochaishoui were closely related. At the end, the new species B. pseudochaishoui Z. Chao sp. nov. was described and illustrated, and a key to the six species was tabulated.

Coupling DNA barcodes and exon-capture to resolve the phylogeny of Turridae (Gastropoda, Conoidea)

Taxon sampling in most phylogenomic studies is often based on known taxa and/or morphospecies, thus ignoring undescribed diversity and/or cryptic lineages. The family Turridae is a group of venomous snails within the hyperdiverse superfamily Conoidea that includes many undescribed and cryptic species. Therefore ‘traditional’ taxon sampling could constitute a strong risk of undersampling or oversampling Turridae lineages. To minimize potential biases, we establish a robust sampling strategy, from species delimitation to phylogenomics. More than 3,000 cox-1 “barcode” sequences were used to propose 201 primary species hypotheses, nearly half of them corresponding to species potentially new to science, including several cryptic species. A 110-taxa exon-capture tree, including species representatives of the diversity uncovered with the cox-1 dataset, was build using up to 4,178 loci. Our results show the polyphyly of the genus Gemmula, that is split into up to 10 separate lineages, of which half would not have been detected if the sampling strategy was based only on described species. Our results strongly suggest that the use of blind, exploratory and intensive barcode sampling is necessary to avoid sampling biases in phylogenomic studies.

Multigene analyses with a broad sampling in Phytophthora and related genera provide evidence for the monophyly of downy mildews

Downy mildews are the most species-rich group of oomycetes, with more than 700 known species. The relationships within the main downy mildew lineages (i.e. the downy mildews with pyriform haustoria, the downy mildews with coloured conidia, and the brassicolous downy mildews) are increasingly well resolved, and 20 well-characterised monophyletic genera have been described. However, their relationships to each other, the various lineages of graminicolous downy mildews, and to the species subsumed in Phytophthora are still unresolved. Recent phylogenomic studies have suggested a polyphyly of the downy mildews, but with a limited taxon sampling within Phytophthora. As taxon sampling is crucial for inferring relationships between large groups, we have conducted a multigene analysis with a set of 72 Phytophthora species and included all known downy mildew lineages. In addition, we performed approximately unbiased (AU) testing as an additional approach to evaluate major nodes. Our analyses resolve the downy mildews as a monophyletic assemblage in all phylogenetic algorithms used. We thus conclude that the evolution of the obligate biotrophy characteristic of downy mildews was a singular event and that all downy mildew pathogens can be traced to a single ancestor.